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In co-annihilation scenarios, the weakly interacting dark matter particle (WIMP) 
is close in mass to another particle that can decay into the WIMP. As a result, the 
other particle does not freeze out before the WIMP in the early universe, and both 
■ contribute to the effective dark matter annihilation cross-section. Since the heavier 

particle does not need to be weakly interacting, the co-annihilation processes are in 
general subject to sizeable radiative corrections. Here this is analyzed for the example 
of neutralino-stop co-annihilation in super symmetry. The leading QCD corrections are 
calculated and it is found that they have a large effect on the effective annihilation 
rS | cross-section, reaching more than 50% in some regions of parameter space. 



In recent years, the existence of dark matter in the universe has been firmly established 
by various different astrophysical experimental methods. From the cosmic microwave back- 
ground and large scale structure, the dark matter density is determined to be f2cDM^ 2 = 
0.1106^o;qq75 [lj, where Ocdm is the ratio of the dark matter energy density to the critical 
density p c = 2Hq / (8tiGn) with the Hubble constant H = hx 100 km/s/Mpc and Newton's 
constant G^. A promising explanation for the nature of dark matter in agreement with direct 
observations and simulations of galaxy formation are weakly interacting massive particles 
(WIMPs). While the Standard Model does not encompass an appropriate particle for this 
purpose, a suitable candidate can arise from many extensions of the Standard Model, such 
as supersymmetry, extended gauge groups or non-anomalous global symmetry groups. 

Due to the requirement that the dark matter particle is weakly interacting, its annihila- 
tion cross-section typically receives only small radiative corrections of 0(%)o This allows 
to make robust predictions of the present dark matter density within a certain model for a 
given choice of parameters. Assuming the standard cosmological model, the evolution of the 
dark matter density from the time of freeze-out from thermal equilibrium until the present 
time can be computed based only on the Boltzmann equation and the thermally averaged 
annihilation cross-section (see for example Ref. [5]). 

Besides predicting a dark matter candidate, most of the Standard Model extensions also 
introduce a large spectrum of additional particles. In supersymmetry, for instance, every 
Standard Model particle has a supersymmetric partner. If the mass of one of these particles 
is close to the WIMP mass, it would not yet be decoupled from thermal equilibrium at 
the freeze-out temperature of the WIMP and thus influences the WIMP annihilation. This 
mechanism is called co-annihilation and effectively lowers the present day relic density [I]. 
It is interesting to note that the co-annihilating particle can have very different quantum 
numbers than the WIMP, as long as it decays into a final state which includes the WIMP. 
In particular, the co-annihilating particle does not need to be weakly interacting. In some 
supersymmetric scenarios, co-annihilation occurs between the lightest neutralino (as the 
WIMP) and the stau or stop, which participate in electromagnetic, and in the latter case 
also strong, interactions. As a result, radiative corrections can be very important in co- 
annihilation processes. 

In this letter, QCD corrections to co-annihilation processes are analyzed. For concrete- 
ness, the specific process of neutralino-stop co-annihilation within the Minimal Supersym- 
metric Standard Model (MSSM) is considered. In many supersymmetry breaking scenarios, 
the light stop t\ is predicted to be the lightest squark state due to large running and mixing 
effects, so that it is not improbable that the t\ can get close in mass to the lightest neutralino 
x\- This scenario is in particular motivated by electroweak baryogenesis, see e.g. Ref. [5]. If 
the mass difference — m^o is small, and all other superpartners are significantly heavier 
than the t\, there are effectively three contributing processes to the evolution of the dark 
matter density, 

xlxl^x, xSti-X, ti*t°-X, (1) 

: An exception are annihilation processes that are suppressed at tree-level, but this suppression is lifted 
at the loop level, see e.g. Ref [2]. 
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Figure 1: Some loop diagrams for the 0(a s ) corrections to (a) X\ti ~ * W + b and (b) 
X% -> tg. 



where X stands for some Standard Model particles. The relic density within a stop-neutralino 
co-annihilation scenario has been first calculated in Ref . [6] , using tree-level formulae for the 
these three annihilation cross-sections. For a typical parameter point with a predominantly 
bino Xi °f H8 GeV and a predominantly R-chiral stop of 138 GeV, the dark matter density is 
predicted to be Vth 2 = 0.112, using DarkSUSY 4.1 [7j. Here XiXi annihilation contributes 
only about 5% to the total averaged annihilation cross-section, while x^ti an d t\t\ annihi- 
lation contribute 85% and 10%, respectively. In the following, QCD radiative corrections to 
the last two processes will be studied in detail. 

X°ti annihilation: The annihilation of x§t\ into Standard Model particles receives con- 
tributions from several channels. For the aforementioned parameter point with Xi ~ B, 
m^o = 118 GeV and t\ ~ t R , = 138 GeV, one finds 

Xiti — > W + b 49.5% contribution to the thermally averaged Xi^i cross-section, (2) 

Xih - tg 47.5%, (3) 

X% - t 7 1-0%, (4) 

Xih -> tZ 2.0%, (5) 

while all other kinematically allowed final states are negligible. Since the first two processes 
are by far dominant, this work concentrates on the radiative corrections to Xi^i ~~ *■ W + b and 
X\ti — > tg. Typical diagrams for the 0(a s ) virtual corrections are shown in Fig. [U Since 
the tree-level processes include quark-squark-neutralino vertices, the QCD corrections with 
gluon exchange cannot be separated from the SUSY-QCD corrections with gluino exchange. 
Therefore the complete SUSY-QCD corrections of order 0(a s ) will be considered. 

The technical methods for the calculation are well established. The virtual loop diagrams 
have been generated with FeynArts 3.0 [8], while the algebraic reduction to standard 
matrix elements was performed with the help of FeynCalc 2.2 [9]. In the loop integrals 
UV divergencies occur, which are canceled through the renormalization of mass, mixing and 
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Figure 2: One diagram in the real gluon corrections to 
Xih — > W + b where the intermediate top-quark propa- 
gator can become on-shell, indicated by a thick line. 

coupling parameters. The masses have been renormalized according to the on-shell scheme, 
while for the strong coupling constant the MS scheme with six-flavor SUSY-QCD running has 
been chosen. Instead of introducing a counterterm for the stop and sbottom mixing angles 
[TO] , the squark sector has been renormalized by using matrix- valued mass counterterms [UJ . 
The two approaches lead to different definitions of the renormalized mixing angles, but are 
equivalent (up to higher orders) when relating physical observables. 

The loop diagrams with gluon exchange also lead to the IR divergencies, that have been 
regularized by a gluon mass. They need to be combined with the contributions with real 
gluon emission in order to arrive at an IR finite complete result. Here a word of caution is 
in order for the real corrections to the process x\t\ ~~ > W + b. In the diagram shown in Fig. [2], 
the top propagator that is shown as a thick line can be become on-shell. In this kinematical 
region, however, it would be a contribution to the process x?^i — * tg with the subsequent 
decay t — > W + b. Therefore the terms with a resonant top propagator need to be subtracted 
from the matrix element of Xiti W + b + g. On the other hand, the interference terms 
between an amplitude with a resonant top propagator and a non-resonant amplitude are 
retained in this channel. 

For the numerical evaluation, the loop integrals are computed with the package Loop- 
Tools [12J, while the phase space for the real corrections is integrated numerically using 
Monte Carlo methods. The real radiation phase space has been mapped onto the integration 
variables in an optimized way so as to improve the result in the regions where the gluon(s) 
become(s) soft or collinear (the collinear case only occurs for the process Xiti tg + g). 
After combining virtual and real corrections, it has been checked that the soft and collinear 
singularities drop out of the total result, so that it does not depend on the gluon mass 
regulator within phase space integration errors. 

In order to simplify the presentation of the numerical results, one can observe that the 
dynamics of co-annihilation are mainly governed by the particle masses, whereas mixing 
effects play a minor role. Thus for simplicity, here it is assumed that the lightest neutralino 
Xi is a pure bino B, while the light stop t\ is purely the partner of the right-handed storJE 
The size of the SUSY-QCD corrections in this scenario is illustrated in Fig. [31 As can be seen 
from the figure, the corrections to Xiti — > W + b amount to only a few per-cent, while the 
process Xiti "~ > tg receives large positive corrections that can reach roughly 50%. Note that 
the size of the corrections is much larger than the scale uncertainty of the Born cross-section, 
see Fig. [U For the next-to- leading order result, the scale uncertainty is only slightly reduced 
with respect to the leading-order cross-section. 

2 In addition, a light stop with a sizeable left-chiral component is in conflict with precision measurements 
of the Zbb vertex. 
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Figure 3: Relative size of the SUSY-QCD corrections to Xih — > W + b (left) and X\ti — > tg 
(right), normalized to the Born cross-sections, respectively. It is assumed that Xi — B and 
t\ — tji. The other relevant supersymmetric parameters are fixed to rriq = 4 TeV for q^t\, 
rrig = 2 TeV, and the QCD scale is set to half the center-of-mass energy, fi Q = y/s/2. 



annihilation: As mentioned above, for small stop masses t\t\ annihilation typically 
contributes much less to the total thermally averaged co-annihilation cross-section then xlti 
annihilation. Nevertheless, there are potentially very large Coulombic gluon corrections that 
can be important for the computation of the relic density. These corrections arise from 
the exchange of gluons between the stop and anti-stop. Since during the phase of freeze- 
out, the stops and anti-stops are slowly moving (E kin j « ^fi. eeze _ out <C m ii)i long-range 
gluon exchange effects lead to a strong enhancement of the tit* annihilation cross-section. 
A similar enhancement effect due to QED corrections in Xi'Xi co-annihilation in the focus 
point region has been studied in Ref. [T3] . 

The leading contribution shown in Fig. [5] (a) leads to a correction 

_ 2 " s7r (R\ 

°coul — — 7, °Born, \S>) 

6V 

where v is the relative velocity between the t\ and t*. The divergence for v — >• is the 
well-known Coulomb singularity. In the case of dark matter annihilation, the Coulomb 
singularity is naturally cut off by the finite temperature in the early universe. Nevertheless, 
v is typically of the same order as a s , so that higher-order effects need to be included. The 
exchange of n gluons, as in Fig. 0(b), generates a correction factor oc (a s /v) n . These effects 
can be systematically resummed within non-relativistic QCD (NRQCD), which is an effective 
theory of full QCD valid at low velocities. Within NRQCD, the problem is described by the 
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Figure 4: Dependence of the cross- 
section a for X\ti —>■ tg on the SUSY- 
QCD scale /i at Born (LO) and one-loop 
(NLO) order. The MSSM parameters 
= B, tt = t R , m^o = 118 GeV, 
138 GeV, m~ q = 4 TeV for q ^ t x , 
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Figure 5: Diagrams for the 
(a ) leading and (b ) resummed 
Coulombic gluon corrections 
to t\t\ annihilation. 



Schrodinger equation with the Coulombic QCD potential V(r) [T 
A 



m 



+ V(r) 



ti 



V(r) = (E + iT)V(r), V(r) 



(7) 



where C-p is the Casimir operator of the SU(3) fundamental representation, and E and V are 
the energy and decay width of the (tit*) system. The correction factor to the cross-section 
is then given by 

0-coul,re S um = |^(q)|2 (g) 



0"Bo 



The QCD Coulomb potential receives important higher order corrections from gluon self- 
interaction and fermion loop effects [T5], [16] . They can be calculated most easily in momen- 
tum space and are given by 



V(q 2 ) = -C F 



47ra s 
q 2 



1 + g (y C A - 2 ^T F n f + ft log(/, 2 /g 2 ) ) + 0(al) 



(9) 



For a typical QCD scale choice /i 



~ m 



ti- 



the logarithm log(/i 2 /g 2 ) = log(l/t>) can become 



large, and it is advantageous to resum contributions of order 0(a s log n v) in a similar way 
as the 0((a s /v) n ) terms. This can be achieved by more elaborate effective theory frame- 
works, such as velocity non-relativistic QCD (vNRQCD) [17] or potential non-relativistic 
QCD (pNRQCD) [TS]. In this work, however, only a simple first estimate of the threshold 
corrections shall be obtained, for which the framework of NRQCD as in eq. ([7J is sufficient. 

The leading contribution to the tit* annihilation cross-section arises from the lowest 
S-wave state, i.e. the (tit*) is state. In eq. ([TJ the decay width T of the stopponium state 
(iit*)is needs to be included. Taking into account the two largest decay channels, (tit*)is 
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Figure 6: Cross-section for t\t\ an- 
nihilation including Coulombic QCD 
corrections near threshold relative to 
the Born cross-section, function 
of the velocity v — (1 — 4m? /s) 1 / 2 . 
The QCD scale has been fixed to fi = 
122 GeV, using five- flavor running of 
a s in normal QCD (not SUSY-QCD). 



gg, W + W , the partial widths at Born level read 
rmths - 99} = 243 l \ mitDis - WW] = 



m- 4m~ 
ti ti 



(10) 



where the second formula is valid in the limits M^ -C 4m^ <C M^ and tan /3 ^> 1. yt is the 
top Yukawa coupling. For mi x = 138 GeV, the total width is T w 5 MeV, and thus almost 
negligible. The stop t\ also has an intrinsic width, with for small mass differences m^o — 
is however smaller than 1 MeV [19] and therefore can be safely ignored. 

For the numerical evaluation, the Schrodinger equation with the QCD potential at two- 
loop order [16] is solved numerically. The result is shown in Fig. [6], which illustrates that 
for v < 0.4 the Coulombic correction can be larger than 100% of the tree-level cross-section. 

For very low values of v, the results of this calculation are not reliable, since the effect of tit^ 
bounds states needs to be taken into account properly. However, as already mentioned above, 
the formation of bound states is inhibited by the non-zero temperature during freeze-out, so 
that the region m^v Tf reez e-out is not relevant for stop annihilation. In other words, in the 
thermally averaged integration over the annihilation cross-section, the region around v = 
has vanishing integration measure, so that no special treatment for the Coulomb singularity 
is necessary. 



Effect on relic density: In order to study the effect of the radiative corrections to 
and tit\ annihilation, as presented in the previous sections, they have been implemented into 
DarkSUSY 4.1 [7|. The effect of the corrections for different neutralino and stop masses 
is shown in Fig. [71 In the left panel only the SUSY-QCD corrections to annihilation 
are incorporated, which can change the predicted relic density by up to about 15%. This 
is comparable with the current uncertainty of ^cdm from astrophysical observations [I]. In 
the right panel, in addition the Coulombic QCD corrections to t\t\ annihilation have been 
taken into account. It can be seen that for small ti-x? mass differences, the predicted value 
of ^cdm can shift by more than 50% due to these effects. Even in the region which is an 
agreement with the current value f2cDM^ 2 = O.llOG^ooyg, the loop correction effects are 
larger than the current error from Ref. pQ. This is indicated by the blue band in the figure, 
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Figure 7: Relative shift of the predicted dark matter density Qcdm due to radiative cor- 
rections. Left: Only SUSY-QCD corrections to Xiti annihilation. Right: SUSY-QCD cor- 
rections to x^ti annihilation and Coulombic QCD to t\t\ annihilation. The blue/dark and 
red/light bands indicate the parameter region in agreement with the observed value of ^cdm 
without and with the radiative corrections, respectively. The parameters have been chosen 
as in Fig. 03 



which is the astrophysically allowed region for tree-level annihilation cross-sections, and the 
red band, which is the same after including the radiative corrections. 

In summary, the leading 0(a s ) corrections to neutralino-stop co-annihilation have been 
calculated. The process is typical for any model and scenario, where a WIMP dark matter 
particle annihilates together with a strongly interacting particle that can decay into the 
WIMP. It was found that the Xit\ annihilation process receives sizeable radiative corrections 
that modify the predicted dark matter relic density by 5-15%. On the other hand, the t\t\ 
annihilation process receives very large corrections, that are associated with the Coulomb 
singularity near threshold. The effect of these Coulombic QCD corrections can have an 
impact on the computed relic density by more than 50% for small stop-neutralino mass 
differences. In a previous work by Hisano et al. [13], similarly large effects have been found 
as a result of QED threshold corrections in neutralino-chargino co-annihilation. 

Since the leading radiative corrections presented in this paper turn out to be large, 
the results are expected to still have a sizeable theoretical error. For a reliable prediction 
of the dark matter relic density, the theoretical calculations need to be refined beyond the 
techniques employed in this work, including the resummation of large logarithms. However, it 
is interesting to note that the effect of the large QCD threshold corrections to tit* annihilation 
can also be derived from a measurement of the pair production cross-section for e + e~ — > t\t\ 
at a future linear collider, see e.g. Ref. [20] . 
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